REPORT DOCUMENTATION PAGE 


Form Approved OMB NO. 0704-0188 


The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, 
searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments 

regarding this burden estimate or any other aspect of this collection of information, including suggesstions for reducing this burden, to Washington 

Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA, 22202-4302. 
Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any oenalty for failing to comply with a collection of 
information if it does not display a currently valid OMB control number. 


PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 


1. REPORT DATE (DD-MM-YYYY) 

23-12-2012 

2. REPORT TYPE 

Final Report 

3. DATES COVERED (From - To) 

1-May-2009 - 30-Apr-2012 

4. TITLE AND SUBTITLE 

Final Report 

5a. CONTRACT NUMBER 

W91 INF-09-1-0111 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

611102 

6. AUTHORS 

Richard A. Regueiro 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

- 1 - 


7. PERFORMING ORGANIZATION NAMES AND ADDRE: 

University of Colorado - Boulder 
The Regents of the University of Colorado 
Office of Contracts and Grants 

Boulder, CO _ 80309 -0572 

9. SPONSORING/MONITORING AGENCY NAME(S) AND 
ADDRESS(ES) 

U.S. Army Research Office 
P.O.Box 12211 

Research Triangle Park, NC 27709-2211 

12. DISTRIBUTION AVAILIBILITY STATEMENT 
Approved for Public Release; Distribution Unlimited 

13. SUPPLEMENTARY NOTES 
The views, opinions and/or findings contained in this report are those of the author(s) and should not contrued as an official Department 
of the Army position, policy or decision, unless so designated by other documentation. 

14. ABSTRACT 

Using solely grain-scale physics-based simulation methods, it is too computationally intensive to account for both 
(I) global initial boundary value problem (IB VP) conditions, and (II) grain-scale material behavior, to understand 
fundamentally the mechanics of dynamic failure in bound particulate materials. The objective of the proposed 
research is to achieve this understanding by accounting simultaneously for grain-scale physics and macro-scale 
continuum IB VP conditions. To achieve the proposed objective, a concurrent computational multi-scale modeling 

15. SUBJECT TERMS 

finite strain, micromorphic, three-dimensional mixed finite element implementation, DE-FE overlap coupling 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 

ABSTRACT 

UU 

15. NUMBER 

OF PAGES 

19a. NAME OF RESPONSIBLE PERSON 

Richard Regueiro 

a. REPORT 

UU 

b. ABSTRACT 

UU 

c. THIS PAGE 

UU 

19b. TELEPHONE NUMBER 

303-492-8026 


Standard Form 298 (Rev 8/98) 


8. PERFORMING ORGANIZATION REPORT 
NUMBER 


10. SPONSOR/MONITOR'S ACRONYM) S) 
ARO 

11. SPONSOR/MONITOR'S REPORT 
NUMBER(S) 

56031-EG.7 


Prescribed by ANSI Std. Z39.18 


































Final Report 


Report Title 


ABSTRACT 

Using solely grain-scale physics-based simulation methods, it is too computationally intensive to account for both (I) global initial boundary 
value problem (IB VP) conditions, and (II) grain-scale material behavior, to understand fundamentally the mechanics of dynamic failure in 
bound particulate materials. The objective of the proposed research is to achieve this understanding by accounting simultaneously for 
grain-scale physics and macro-scale continuum IBVP conditions. To achieve the proposed objective, a concurrent computational 
multi-scale modeling approach will be developed that involves the following 3 features: (1) coupling regions of micromorphic continuum 
finite element to an 'open window' on the particulate micro-structure where localized deformation nucleates and an interface with a 
deformable solid body could exist; (2) converting to discrete element fragmentation modeling in micro-structural regions; and (3) adapting 
numerically grain-scale resolution over the material domain. The desired result is to enable a more complete understanding of the role of 
grain-scale physics on the thermo-mechanical properties and performance of heterogeneous bound particulate materials of interest to the 
Army. The contribution will be demonstrated by working with Army researchers in the Impact Physics Branch during onsite visits and 
through an open-source research code Tahoe where the concurrent multi-scale model will be implemented. 

Enter List of papers submitted or published that acknowledge ARO support from the start of 
the project to the date of this printing. List the papers, including journal references, in the 
following categories: 

(a) Papers published in peer-reviewed journals (N/A for none) 


Received Paper 


08/29/2011 2.00 R. A. Regueiro, V. Isbuga. Length scale effects in finite strain micromorphic linear isotropic elasticity: finite 

element analysis of three-dimensional cubical microindentation, 

Proceedings of the Institution of Mechanical Engineers, Part N: Journal of Nanoengineering and 
Nanosystems, (07 2011): 0. doi: 10.1177/1740349911412850 

08/29/2011 3.00 Volkan Isbuga, Richard A. Regueiro. Three-dimensional finite element analysis of finite deformation 

micromorphic linear isotropic elasticity, 

International Journal of Engineering Science, (05 2011): 0. doi: 10.1016/j.ijengsci.2011.04.006 
TOTAL: 2 


Number of Papers published in peer-reviewed journals: 

(b) Papers published in non-peer-reviewed journals (N/A for none) 


Received Paper 


TOTAL: 





Number of Papers published in non peer-reviewed journals: 


(c) Presentations 

October 2012, Concurrent multiscale computational modeling of particulate materials: (i) 3D ellipsoidal DEM of X-ray CT sand 
assemblies, (ii) overlap coupling of DEM and micromorphic FEM, SES 2012, GeorgiaTech, Atlanta. 

September 2012, Total Lagrangian Nonlinear Finite Element Analysis of Finite Strain Micromorphic Pressure-Sensitive Elastoplasticity, 
ECCOMAS 2012, University of Vienna, Austria. 

Number of Presentations: 2.00 

Non Peer-Reviewed Conference Proceeding publications (other than abstracts): 


Received Paper 


12/23/2012 5.00 Volkan Isbuga, Richard Regueiro. FINITE ELEMENT ANALYSIS OF FINITE STRAIN MICROMORPHIC 

ISOTROPIC ELASTICITYAND DRUCKER-PRAGER PLASTICITY, 

2012 ASCE Joint Conference of the Engineering Mechanics Institute. 2012/06/18 02:00:00,. : , 

TOTAL: 1 


Number of Non Peer-Reviewed Conference Proceeding publications (other than abstracts): 


Peer-Reviewed Conference Proceeding publications (other than abstracts): 


Received Paper 


08/29/2011 i .00 Richard A. Regueiro, Beichuan Yan. Coupling discrete elements and micropolar continuum through an 
overlapping region, 

GeoFrontiers 2011.2011/03/14 02:00:00,. : , 

TOTAL: 1 


Number of Peer-Reviewed Conference Proceeding publications (other than abstracts): 


(d) Manuscripts 







Received 


Paper 


TOTAL: 


Number of Manuscripts: 


Books 


Received Paper 


12/23/2012 6.00 Richard Regueiro, Beichuan yan. Computational homogenization and partial overlap coupling between 

micropolar elastic continuum finite elements and elastic spherical discrete elements in one dimension, 
Boca Raton, FL: Pan Stanford Publishing Co., (02 2013) 

TOTAL: 1 


Patents Submitted 


Patents Awarded 


Awards 


Graduate Students 


NAME 

PERCENT SUPPORTED Discipline 

Volkan Isbuga 

0.50 

FTE Equivalent: 

0.50 

Total Number: 

1 


Names of Post Doctorates 


NAME 

PERCENT SUPPORTED 

Volkan Isbuga 

0.10 

FTE Equivalent: 

0.10 

Total Number: 

1 


Names of Faculty Supported 















NAME 

PERCENT SUPPORTED National Academy Member 

Richard Regueiro 

0.08 

FTE Equivalent: 

0.08 

Total Number: 

1 

Names of Under Graduate students supported 

NAME 

PERCENT SUPPORTED 

FTE Equivalent: 


Total Number: 



Student Metrics 

This section only applies to graduating undergraduates supported by this agreement in this reporting period 

The number of undergraduates funded by this agreement who graduated during this period:. 0.00 

The number of undergraduates funded by this agreement who graduated during this period with a degree in 

science, mathematics, engineering, or technology fields:. 0-00 

The number of undergraduates funded by your agreement who graduated during this period and will continue 

to pursue a graduate or Ph.D. degree in science, mathematics, engineering, or technology fields:. 0-00 

Number of graduating undergraduates who achieved a 3.5 GPA to 4.0 (4.0 max scale): 0.00 

Number of graduating undergraduates funded by a DoD funded Center of Excellence grant for 

Education, Research and Engineering:. o.OO 

The number of undergraduates funded by your agreement who graduated during this period and intend to 

work for the Department of Defense. 0.00 

The number of undergraduates funded by your agreement who graduated during this period and will receive 

scholarships or fellowships for further studies in science, mathematics, engineering or technology fields:. 0.00 

Names of Personnel receiving masters degrees 

NAME 

Total Number: 

Names of personnel receiving PHDs 

NAME 
Volkan Isbuga 

Total Number: 1 

Names of other research staff 

NAME PERCENT SUPPORTED 

FTE Equivalent: 

Total Number: 

Sub Contractors (DD882) 



























Inventions (DD882) 


See attachment 


Scientific Progress 


Technology Transfer 



Final Report for W911NF-09-1-0111: period 05/01/2009 to 09/30/2012 


1 


Scientific Progress and Accomplishments 

contract /grant number: W911NF-09-1-0111 
ARO proposal number: 56031-EG 
December 31, 2012 
Richard A. Regueiro 
Associate Professor 

Department of Civil, Environmental, and Architectural Engineering 
University of Colorado, Boulder 

The outline of this section of the final report is as follows: 

1. overview of micromorphic finite element (FE) results presented in various published 
journal articles [Isbuga and Regueiro, 2011, Regueiro and Isbuga, 2011, Isbuga and 
Regueiro, 2012] and PhD thesis [Isbuga, 2012] funded by this research grant; 

2. overview of discrete element and micropolar continuum finite element coupling in one 
dimension, and numerical results [Regueiro and Yan, 2013]; 

3. formulation of overlap coupling for grain-FE scale and micromorphic continuum FE 
scale (two methods); 

4. overview of adaptivity approach; 

5. summary of future research 

1 Overview of micromorphic finite element (FE) re¬ 
sults 

This section summarizes briefly the micromorphic finite element (FE) results presented in 
various published journal articles [Isbuga and Regueiro, 2011, Regueiro and Isbuga, 2011, 
Isbuga and Regueiro, 2012] and PhD thesis [Isbuga, 2012] funded by this research grant. 
More details can be found in the cited articles and thesis, provided as attachments. 

The papers Isbuga and Regueiro [2011], Regueiro and Isbuga [2011] focussed on the finite 
strain, materially linear isotropic clastic, three-dimensional finite element implementation 
and numerical results using Tahoe. We showed that the implementation is convergent with 
respect to spatial discretization (i.e., mesh size), even in the presence of a boundary layer. 
Upon increasing elastic length scale, the boundary layer was found to become less sharp. 
For pressure-sensitive Drucker-Prager elasto-plasticity [Isbuga and Regueiro, 2012, Isbuga, 
2012], it was found that choice of micromorphic elastic parameters in conjunction with plastic 
strength parameters will influence the stress paths and, in turn, yielding. The boundary 
conditions on the micro-displacement tensor had a significant influence on the results. A 
plastic length scale term has not yet been implemented, and it is anticipated that this will 
influence the numerical results where there is a gradient of an internal state variable (ISV) 
in a region of large deformation. 
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2 Overview of discrete element and micropolar contin¬ 
uum finite element coupling in one dimension, and 
numerical results 

The ultimate problem of interest is shown in Fig.l. The discrete element (DE) and/or finite 
element (FE) representation of the particulate micro-structure is intentionally not shown 
in order not to clutter the drawing of the micro-structure. The grains (binder matrix not 
shown) of the micro-structure are ‘meshed’ using DEs and/or FEs with cohesive surface el¬ 
ements (CSEs). The open circles denote continuum FE nodes that have prescribed degrees 
of freedom (dofs) D based on the underlying grain-scale response, while the solid circles 
denote continuum FE nodes that have free dofs D governed by the micromorphic continuum 
model. We intentionally leave an ‘open window’ (i.e., Direct Numerical Simulation (DNS) 
region) on the particulate micro-structural mesh in order to model dynamic failure. If the 
continuum mesh overlays the whole particulate micro-structural region, as in Klein and Zim¬ 
merman [2006] for atomistic-continuum coupling, then the continuum FEs would eventually 
become too deformed by following the micro-structural motion during fragmentation. The 
blue-dashed box at the bottom-center of the illustration is a micromorphic continuum FE 
region that can be converted to a DNS region for adaptive high-fidelity material modeling 
as the projectile penetrates the target. 
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Figure 1. 2D illustration of concurrent computational multi-scale modeling approach in the contact 
interface region between a bound particulate material (e.g., ceramic target) and deformable solid 
body (e.g., refractory metal projectile). 
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In order to simplify the problem, we consider the compression and shear of the interface 
of DNS region with the micromorphic continuum FE region as a ID overlap problem between 
glued, nonlinear elastic DEs, and micropolar elastic continuum FE [Regueiro and Yan, 2013]. 
Details may be found in Regueiro and Yan [2013]. In summary, we learned that the coupling 
methodology generally works for quasi-statics between DEs and FEs, but the ID problem 
is too “discrete” to yield the results we expect, i.e., no influence of overlap coupling region. 
This is because shear and rotational degrees of freedom are coupled in the DE and FE 
formulations, and thus upon shearing the end of the DE string of particles, a ratcheting 
rotational mode is propagated through the DE string and into the FE continuum domain 
where it is smoothed out. We are working on extending this coupling method to include 
inertia terms (dynamics), and also to couple grain-FE to micromorphic FE for 2D and 3D 
problems. 
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3 Formulation of overlap coupling for grain-FE scale 
and micromorphic continuum FE scale 

Two coupling approaches in Sections 3.1 and 3.2 were considered in this research. Cou¬ 
pling Approach 1 is implemented with numerical results discussed in Section 2. Coupling 
Approach 1 is based on the “bridging-scale decomposition” in Wagner and Liu [2003], Kad- 
owaki and Liu [2004], Klein and Zimmerman [2006], and Coupling Approach 2 is based on 
the “bridging domain method” in Xiao and Bclytschko [2004]. Coupling Approach 1 couples 
the DNS and micromorphic FE dofs directly through interpolation functions in the overlap 
region, and partitions the energy through scaling factors within each term of the energy 
(kinetic and potential energies). Coupling Approach 2 partitions the energies through a 
scalar distance function ce(JC) (describing the width of the overlapping domain), and uses 
Lagrange multipliers to constrain the DNS displacements to the micromorphic continuum 
FE displacements in the overlapping domain. 

3.1 Coupling Approach 1: 

It is assumed that the micromorphic contimium-FE mesh covers the domain of the problem 
in which the bound particulate mechanics is not dominant, whereas in regions of significant 
grain-matrix debonding or intra-granular cracking leading to a macro-crack, a grain-scale 
mechanics representation is used (e.g., grain-FE, grain-DE-FE, grain peridynamics [Silling, 
2000]) in Fig.2. 

3.1.1 3D Kinematics 

Here, a summary of the kinematics of the coupled regions is given for general 3-D kinematics, 
following the illustration shown in Fig.2. It is assumed that the micromorphic continuum 
finite element (FE) mesh covers the domain of the problem in which the material is not 
fragmenting, whereas in regions of fracture and fragmentation, a DNS (e.g., DE/FE/CSE, 
or peridynamics [Silling, 2000]) representation can be used. In Fig.2, discrete domains are 
defined, where the yellow background denotes the FE overlap region B h with underlying 
ghost DNS nodes, brown the micromorphic FE continuum region B h with no underlying 
DNS, and white background the free DNS region B h UB DNS . In summary, the finite element 
domain B h is the union of pure micromorphic continuum FE domain B h , overlapping FE 
domain with underlying ghost DNS nodes B h , and overlapping FE domain with underlying 
free DNS nodes B h , such that B h = B h UB h UB h . The pure DNS domain with no overlapping 
FE domain (i.e., the ‘open-window’) is indicated by B DNS . The goal is to have the overlap 
region B h U B h as close to the region of interest (e.g., penetrator skin) as to minimize the 
number of DNS nodes, and, thus, computational effort. Following some of the same notation 
presented in Klein and Zimmerman [2006], we define a generalized dof vector Q for DNS 
displacements in the system as 

Q = [q a , Qp, ■ ■ ■ , <Z 7 ] r , a,P,...,'yeA (1) 

where q a is the displacement vector of DNS node a, and A is the set of all DNS nodes. 
Likewise, the micromorphic finite element nodal displacements and micro-displacements are 
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Figure 2. Two-dimensional illustration of the Coupling Approach 1 between grain-FE DNS and 
micromorphic FE continuum regions. 


written as 

D = [d a ,d b , ...,d c ,0 d ,0 e ,...,0 / ] T (2) 

a,b,... ,c E J\ , d, e,..., f e M. 

where d a is the displacement vector of node a , <p d is the micro-displacement tensor in vector 
form of node d, N is the set of all nodes, and A4 is the set of micromorphic finite element 
nodes with micro-displacement tensor degrees of freedom, where A4 C JV. In order to satisfy 
the boundary conditions for both regions, the motion of the DNS nodes in the overlap region 
(referred to as “ghost DNS nodes,” cf. Fig.2) is prescribed by the micromorphic continuum 
displacement and micro-displacement fields, and written as 

Q = [q a , Q/}, ■ ■ ■, q-y] T , a,/3,...,7 e A, Ae B h (3) 

while the unprescribed (or free) DNS nodal displacements are 

Q = \qs,q„----%} T , Ac& h uB DNS ( 4 ) 
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where A U A = A and A D A = 0. Likewise, the displacements and micro-displacement 
tensor components of nodes overlaying the DNS region are prescribed by the DNS motion 
and written as 


D = [d a ,d b , ...,d c ,(j) d ,(t> e ,...,(j) f ] T (5) 

a, b ,..., c e Af , d,e ,..., / e M 
Af,M e B h u B h 

while the nnprescribed (or free) nodal displacements and micro-displacements are 

D = [d m , d n ,..., d s , (fit, (j) u , ..., <fi v \ (6) 

m,n,..., s E Af , t, u ,..., v E A4 
Af,Me B h u B h 

where J\f U J\f = J\f, Af C\ J\f = (ft, M U M = M, and A4 ft A4 = 0. Referring to Fig.2, 
the prescribed DNS nodal dofs Q can be viewed as boundary constraints on the free DNS 
region, and likewise the prescribed micromorphic finite element nodal displacements and 
micro-displacements D can be viewed as boundary constraints on the finite element mesh 
in the overlap region. 

In general, the displacement vector of a DNS node a can be represented by the fi¬ 
nite element interpolation of the micromorphic continuum macro-displacement field u h and 
micro-displacement tensor held <& h evaluated at the DNS node in the reference configuration 
X^, such that 

= £KPdKW , = ENPCMW aeA (7) 

aeAf b£M 

where N/ are the shape functions associated with the micromorphic continuum displacement 
held u h , N b are the shape functions associated with the micromorphic continuum micro¬ 
displacement held <& ', and h indicates the characteristic length of the micromorphic FE 
mesh. Recall that N/ and Nf have compact support and, thus, are only evaluated for DNS 
nodes that lie within a micromorphic FE containing nodes a and b in its domain. 

We now derive the micro-continuum displacement (u') h of a micro-element (Fig.3) using 
the interpolated macro-displacement u h and micro-displacement tensor <f> h in Eq.(7). The 
position of DNS node a in the current configuration is 

(x') h (X h a , E h a , t ) = x h (X h Q , t ) + e\X h a , S„, t) (8) 

where is the reference relative position vector of DNS node a with respect to reference 
continuum position X h a (e.g., a Gauss point) in Bq, and $, h is the current relative position 
vector of DNS node a with respect to the current continuum position x h in B h . Introducing 
the relations, 


* {X a, *) 
(x h a ,s h a ,t) 


X h a + u h (X h a ,t) 

X h (X h a ,t)-~ h a = (l + <f> h (X h a ,t))-Z h a 


and substituting into Eq.(8), we have the current position of DNS node a as 


(9) 

( 10 ) 
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Figure 3. Multiplicative decomposition of deformation gradient F and micro-deformation tensor 
X • Geometrical points (“macro-elements”) with centroids C, C , and c live in their respective 
configurations: P £ £>o, P £ B, and p £ B. Material points (“micro-elements”) with centroids 
C', C\ and d. Differential vectors and deformable directors are mapped accordingly: dx = FdX, 
dx = F e dX , dX = F p dX , £ = * 3 , £ = x e S, and S = X P S. 

(xTiK, SS, *) = X a + u\ x l t) + (1 + * h (X h a: t)) ■ S h a (11) 

If we rewrite {x') h in terms of its reference position ( X') h , and we introduce a micro¬ 
continuum displacement vector ( u') h , we have 

{x') h \X h al t ) = (. X')\X h a , Sq) + {u')\Xl S h a , t) (12) 

where ( X') h (X h a , E h a ) = X h a + E h a , then we have 

(u') h (X h a , Sai t) = C x')\X h a , Si, t) - Xj - (13) 

Finally, we have the micro-continuum displacement vector ( u') h associated with DNS node 
a, such that 

q a (t) = (u') h (X h a ,~ h a ,t) = u h (X h a ,t ) + * h (X h a ,t) ■ S h a (14) 

Then, for all ghost DNS nodes (cf. Fig.2), the interpolations can be written as 

Q = Xq D ■ D + ■ D (15) 

where Nq D and Ngg are shape function matrices containing individual nodal shape func¬ 
tions N/ and N®, but for now these matrices will be kept general to increase our flexibility 
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in choosing interpolation/projection functions (such as those used in meshfree methods). 
Overall, the DNS nodal displacements may be written as 

Q 

Q 

where Q' is introduced [Klein and Zimmerman, 2006] as the error (or “fine-scale” [Wagner 
and Liu, 2003]) in the interpolation of the free DNS nodal displacements Q, whose function 
space is not rich enough to represent the true free DNS motion. The shape function matrices 
N are, in general, not square because the number of free DNS dofs are not the same as free 
and prescribed micromorphic FE dofs, and the number of ghost DNS dofs is not the same 
as prescribed and free micromorphic FE dofs. A scalar measure of error e in DNS nodal 
displacements is defined as [Klein and Zimmerman, 2006] 


N 


QD 


N 


QD 


n 6d n. 



D 

+ 

Q' 


D 


0 


(16) 


e = Q' ■ Q' 


(17) 


which may be minimized with respect to prescribed micromorphic continuum FE nodal dofs 
D to solve for D in terms of free DNS dofs and micromorphic continuum FE nodal dofs as 

5 = M BB N Id(Q - N QoD) . M m = N t qB N qB (18) 

This is known as the “discretized L 2 projection” [Klein and Zimmerman, 2006] of the free 
DNS nodal dofs Q and free micromorphic FE nodal dofs D onto the prescribed micromorphic 
FE nodal dofs D. Upon substituting Eq.(18) into Eq.(15), we may write the prescribed 
DNS dofs Q in terms of free DNS dofs Q and micromorphic continuum FE nodal dofs D. 
In summary, these relations are written as 


where 


Q 

— BqqQ + Bq D D 

(19) 

D 

= BgqQ + b Bd d 

(20) 

B QQ 

= Nq B B B q 

(21) 

B QD 

— Nqd + NqqB Bd 

(22) 

B DQ 

= MWivW 

DD QD 

(23) 

■®6o 

= -M^N^Nq D 

DD QD ^ 

(24) 


As shown in Fig.2, for a finite element implementation of this dof coupling, we expect that 
free DNS dofs Q will not fall within the support of free micromorphic continuum FE nodal 


dofs D , such that it 

can be assumed that Nqd = 0 and 



Q = B Bq Q + Bq d D , D = BqqQ 

(25) 

where 

Bqq = Nq B Bqq , Bq d = Nq d 

(26) 


B dq = M bb N t qB , B Bd = 0 

(27) 


The assumption Nqd / 0 would be valid for a meshfree projection of the DNS nodal dofs 
to the micromorphic FE nodal dofs, as in Klein and Zimmerman [2006], where we could 
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imagine that the domain of influence of the meshfree projection could encompass a free 
DNS node; the degree of encompassment would be controlled by the chosen support size of 
the meshfree kernel function. The choice of meshfree projection in Klein and Zimmerman 
[2006] was not necessarily to allow Q be projected to D (and vice versa), but to remove the 
computationally costly calculation of the inverse IVT-h in Eqs. (19) and (20). 

3.1.2 3D Kinetic and Potential Energy Partitioning and Coupling 

For the DNS equations, the kinetic energy is T®, dissipation function F®, and potential 
energy U®, such that 


T q = i QM q Q 

F Q = a Q T Q (28) 

rQ 

U q (Q) = / F 1NT,Q (S)dS 

Jo 

where is the mass matrix, and F INT ’® is an internal force vector for the DNS region. 
The dissipation function F® is written as a linear function of the kinetic energy T® with 
proportionality coefficient czq, which falls within the class of damping called “Rayleigh damp¬ 
ing” (pg.130 of Rayleigh [1945]). For the micromorphic continuum FE equations, T D is the 
kinetic energy, F D the dissipation function, and U D the potential energy, such that 


rj~\D _ 


F d = 


-DM d D 

2 

a D T D 

,D 


U d (D) 


? INT,D 


1 0 


(S)dS 


(29) 


where M D is the mass and micro-inertia matrix, and F 1 ;V/ ' I} is the internal force vec¬ 
tor. We assume the total kinetic and potential energy and dissipation of the coupled DNS- 
micromorphic-continuum system may be written as the sum of the energies 


T(Q,D) 

= r«(Q,Q(Q,D)) + T D (D,D(Q)) 

(30) 

U(Q,D) 

= UQ(Q,Q(Q j D)) + U d (D.D(Q)) 

(31) 

F(QD ) 

= F q (Q.^(Q,D)) + F d (D,D(Q)) 

(32) 


where we have indicated the functional dependence of the prescribed DNS dofs and mi¬ 
cromorphic FE dofs solely upon the free DNS dofs and micromorphic FE dofs Q and D. 
respectively. Lagrange’s equations may then be stated as 


d fdT\ dT OF dU 
dt \dQ ) ~ dQ + dQ + dQ 
d (dT\ dT dF dU 
~dt\dDj ao + ao + ao 


piEXT.Q 

pEXT,D 


(33) 
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which lead to a coupled system of governing equations (linear and first moment of momen¬ 
tum) for the coupled DNS-micromorphic-continuum mechanics. The derivatives are 


dT 

dQ 

dF 

dQ 

dU 

dQ 

dT 

&d 

dF 

dt> 

dU 

dD 


dT Q dT Q 


dT D 


„ dT n 
’ dQ~° 



an 


= 0 


st d 

dD 


B 


DQ 


fST D \ 

\~8d) 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


If the potential energy U is nonlinear with regard to DNS constitutive models and rnicro- 
morphic elasto-plasticity, then Eq.(33) may be integrated in time and linearized for solution 
by the Newton-Raphson method. The benefit of this multiscale method, as pointed out by 
Wagner and Liu [2003], is that time steps may be different for the DNS and micromorphic 
FE solutions. A multiscale time stepping scheme will follow an approach similar to Wagner 
and Liu [2003]. To complete Eq.(33) and identify an approach to energy partitioning, the 
individual derivatives may be written as 


dT Q 

~dQ 

dT D 

dD 

dF 

db 

dF 

~dD 

dU Q 

~dQ 

dU D 

dD 


ary - a 

= M q Q , —= M q Q 

dQ 

dT D ~ a 

= m' d , —= M d D 

dD 

a Q (M«Q + + a D 


a Q (B T dD M^Q]+a D (m d D 


QD 

F int ^(Q) , = F int, ®(Q) 

dQ 

pi NT,Did) , dC = F INT -' D {b) 

OD 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


where superscript Q denotes free DNS dofs and Q ghost DNS dofs, whereas superscript D 
denotes free micromorphic FE dofs and D prescribed micromorphic FE dofs. The energy 
partitioning will be introduced through the definition of these terms below. First, substitute 
Eqs. (40)-(45) into Eq.(33) to arrive at the coupled nonlinear equations in terms of Q and 
D as 




































Final Report for W911NF-09-1-0111: period 05/01/2009 to 09/30/2012 


11 


(m« + B t Qq M$B 3q + Bl Q M B B BQ ) Q 
+ ( B S Q mQb Gd + b I q M B B S d ) D 

+ /qM® + a + a D B Bq M~ B B( fj Q + (iqBqqM®Bq D D 
+F ,nt *(Q) + B t Sq F' nt - B {B S q Q + b Sd d} 

+B t Bq F ,nt - S IB B q Q + B Bd D] = F ext -0 + B t Sq F ext -° (46) 

+ B t Bd M B B B q ) Q 

+ (m d + b| d m«b 6d + b t Bd m b b B d ) D 

+ a Q B q D M B BqqQ + (a 0 B gD M«B 6D + a B M j D 

+B t Sd F ,ntB [B S q Q + Bq d D] + B t Bd F ,nt ' d [B B q Q + B Bd D] 

+F ,nt ’ d (D) = F ext - d + b t Bd F extB (47) 

where an expression in brackets [•] denotes the functional dependence of the nonlinear inter¬ 
nal force vector. Note the projections through the B matrices of the corresponding mass and 
damping matrices, and forcing vectors. First, starting with the mass matrices for the DNS 
nodes, we partition the kinetic energy (and, in turn, the dissipation functions) as follows: 

M® = (1 - f ) A m a , aei, x a e B h (48) 

a 

M q = M dns ’ q + M Q (49) 

m dns,q = , p e A , Xp e B DNS 

P 

M Q = (1 — f) ^ m/p , /3 e A , e B h 

0 

where is the mass matrix of ghost DNS nodes in B h , M DNS ’® the mass matrix of free 

DNS nodes in B DNS , M the mass matrix of free DNS nodes in B h , and f is a weighting 
factor for the kinetic energy in the overlap region B h U B h . For no homogenized micromor- 
phic continuum contribution to the kinetic energy in the overlap region, f = 0, and for full 
micromorphic continuum homogenization of the underlying DNS kinetic energy, f — 1. In 
our case, we will consider the range 0 < f < 1. Given that the proposed multiscale modeling 
framework is to be used in an adaptive fashion in the future, having an overlaying micromor¬ 
phic continuum homogenization of the DNS response is attractive when DNS is converted to 
micromorphic continuum representation, and vice versa (in a statistical manner*). For the 

* statistical, in the sense of generating a DNS representation from a micromorphic continuum one, where 
the underlying DNS region does not exist; converting from DNS to micromorphic continuum representation 
is straightforward given the built-in homogenization that the micromorphic continuum possesses 
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mass and micro-inertia matrices associated with the micromorphic continuum, we partition 
the kinetic energy as follows: 


M d = 


— D - D 

M +M 


(50) 


M = 


M 0 




0 M 

M" = J\(f (m u,e ) + f e m u ’ e ) 


e&B h 


M — (f + f e m #,e ) 


eg B h 


-—D 

M = f 


M 0 


0 M 


M d = 


'—' D — n 

M +M 


M D = 


M 0 
0 M* 


, M U = ^ (m u,e ) AI J\ <m # ’ e > 

eg B h eg B h 


, M“=^ m ue , = A m$,e 

egS h egB^ 1 


(51) 


where IVf is the micromorphic continuum mass and micro-inertia matrix associated with 

~ — D 

prescribed micromorphic FE dofs in £> , M the micromorphic continuum mass and micro- 

— 5 -—D 

inertia matrix associated with free micromorphic FE dofs in B h 1 where M and Af are 

extracted from the total mass and micro-inertia matrix M in B h , with superscript {•)^ > 
denoting the full mass and micro-inertia matrix associated with elements in B h , (•) is a 
homogenization operator, f e is the partitioning coefficient of micromorphic continuum kinetic 
energy associated with element B e C B h . A simple choice is a volume fraction f e = B e,D /B e , 
where B e,D = B e — B e ’®] B/ ,D is the non-overlapping micromorphic continuum part of element 
volume B e C B h , and B e ’® is the overlapped ghost DNS volume in the element (cf. Fig.2). For 
kinetic energy partitioning, a volume fraction that directly relates to mass and micro-inertia 

.— D 

partitioning seems an appropriate choice. M is the homogenized micromorphic continuum 
mass and micro-inertia matrix associated with prescribed micromorphic FE dofs in B h \ where 
iff = 0, there is no micromorphic continuum homogenization in B h (i.e., all kinetic energy is 
due to underlying DNS). is the micromorphic continuum mass and micro-inertia matrix 
associated with free micromorphic FE nodal dofs in the pure micromorphic continuum FE 
domain B h . 

For the potential energy (internal force) and external force partitioning in the DNS 
system, we write 
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jpINT,Q _ 

(1 - q) A fl NT ’ Q 
£ 

(52) 

j?INT,Q _ 

jpINT,DNS,Q _j_ f INT ’Q 

(53) 


F INT,DNS,Q = A f INT,Q xs e B DNS 
5 


F INT,Q = (1 - q) A fs NT ’ Q ,x s eB 



<5 


f EXT,Q 

= (i-g)A ff XT ' Q ,x e eB h 

(54) 

f EXT,Q 

6 

pEXT,DNS,Q f EXT ’Q 

(55) 


f ext,dns,q = A fEXT.Q Xs e b dns 
§ 



F EXT ' Q = (1 - q) A f! XT ’ Q , x 5 e B h 



5 


where F 1NT ’Q i s the internal force vector associated with ghost DNS nodes in B h , jt int Fns,q 

^int,q 

is the internal force vector associated with free DNS nodes in B DNS , F is the internal 

A j 771 ym N 

force vector associated with free DNS nodes in B n , F ,w is the external force vector asso- 
dated with ghost DNS nodes in B n , F ’ is the external force vector associated with 
free DNS nodes in B UNb , F is the external force vector associated with free DNS nodes 
in B h , and q is a weighting factor for the potential energy in the overlap region B h U B h . For 
no homogenized micromorphic continuum contribution to the potential energy in the overlap 
region, q — 0, and for full micromorphic continuum homogenization of the underlying DNS 
potential energy, q — 1. In our case, we will consider the range 0 < q < 1. Note that in Klein 
and Zimmerman [2006], they chose q = 0. Their Cauchy-Born elastic constitutive model acts 
like a homogenization operator on the underlying atomistic response, but instead of keeping 
an overlain Cauchy-Born representation, the potential energy is completely represented by 
the underlying atomistic response, except in the overlap region B h where partitioning occurs. 
Note that x e and x$ denote positions of DNS nodes for calculating internal force vectors in 
Eqs. (52)-(56), whereas x a and xp in Eqs. (48) and (49) denote DNS nodes for calculating 
DNS mass matrices. 

For the potential energy (internal force) partitioning in the micromorphic continuum, we 
write 


f int,d 


~INT,D p int,D 


~INT,D 

F 


f-,INT,u 

F 

F 


-,INT,u 


A 

e&B h 


j?I NT,u,e 


F 


A/ 


INT,3>,e 


e&B h 


(56) 


and 
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jpINT,D 


~ INT,D — INT,D 

F +F 


(57) 


- INT,D 


- INT,u 


• INT,u 


- /TVT,# 


A 0(f 

e&B h 


INT,u,e 


> + q*f INT ’ u ’ e ) 


INT,u,e\ _ 


f 


■ I NT,® 


> = /(s 


u,e\T 


A 

e£B h 

— I NT,u 


<T) dv 

INT,<fr,e\^ _|_ ~e 


^ I NT,D 

F — q 


-1 NT,u 


A 


— INT,® 

F 

^INT,u,e\^ 


INT,$ 


A </ 

e&B h 


INT,®,e\ 


~INT,D 

where is the internal force vector associated with free micromorphic FE dofs in £1 , 

~ INT,D 

E 1 the internal force vector associated with prescribed micromorphic FE dofs in B tl , 

~ I NT,D ~INT,D ~ INT,D 

where -F and F are extracted from the full internal force vector F , with 
superscript denoting the full internal force vector associated with elements in B h , g e 
is the partitioning coefficient of micromorphic continuum potential energy associated with 
element B e C B h , and (•) is a homogenization operator (to be defined later). A simple choice 
is a volume fraction q e = f e . Klein and Zimmerman [2006] considered a more sophisticated 
approach using an atomic bond density function solved to reproduce a minimum potential 
energy state for homogeneous deformation. The analogy here for DNS nodes would be a DNS 
potential energy term (internal force vectors). This will be considered further in future work. 
For now, we consider a volume fraction partitioning through q e , and a simple scaling through 

— INT,D 

coefficients q. F is the homogenized internal force vector associated with prescribed 
micromorphic FE dofs in B h , which has no contribution if q — 0, i.e., underlying DNS internal 
forces provide full contribution in B h F ’ is the internal force vector associated with 
free micromorphic FE dofs in the pure micromorphic continuum domain B h . The external 
force vector is written as 


jEXT,D 


~EXT,D pEXT,D 


EXT,D 


EXT,u 


-.EXT,® 


p EXT,u = Ft + F u + ^ f 

e£B h 


EXT,u,e 

b 


f ext,® = 


A f 

e£B h 


EXT,®,e 

i 


(58) 
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and 


jp EXT,D 


~EXT,D — EXT,D 

F + F 


- EXT,D 


- EXT,u 


- EXT,® 


■ EXT,u 


- EXT,® 


A U(/ 6 

e&B h 


EXT,u,e 


A 

e&B h 


?(// 


i ~e nEXT.u.e 

+ 9 fb 


+ 9 e ft 


— EXT,D 

F — q 


- 


-EXT,® 


- EXT,u 


A </ 

e£B h 


EXT,u,e 

b 


-EXT,® 


A 

e£B h 


(59) 


/ 


EXT,®,e 

l 


~ EXT,D 

where F is the external body force and couple vector associated with free micromorphic 

~EXT,D 

FE nodal dofs in B h , F the external body force and couple vector associated with 

~ EXT,D ~ EXT,D 

prescribed micromorphic FE dofs in B , where F and F are extracted from 

~ EXT,D ^ EXT,D 

F , the total external body force and couple vector calculated in B h . F is the 

homogenized external body force and couple vector associated with prescribed micromorphic 
FE dofs in B h , which has no contribution if q — 0, i.e., underlying DNS body forces provide 
full contribution in B h . F is the external force and couple vector associated with free 
micromorphic FE dofs in the pure micromorphic continuum FE domain B h . 


3.2 Coupling approach 2: 

In this coupling approach, it is also assumed that the micromorphic continuum-FE mesh 
covers the domain of the problem in which the bound particulate mechanics is not dominant, 
whereas in regions of significant grain-matrix debonding or intra-granular cracking leading 
to a macro-crack, a grain-scale mechanics representation is used in Fig.4. 

In Fig.4, discrete domains are defined, where the yellow background denotes the DNS- 
FE overlap region J3 DNS ~ FE ; with DNS region denoted by B DNS terminating at dB DNS , and 
micromorphic continuum FE region denoted by B h E terminating at dB b E . 

3.2.1 3D Kinematics 

The notation here for Coupling Approach 2 is similar to that used for Coupling Approach 
1. We define a generalized dof vector Q for all DNS nodal displacements in the system as 

Q = [q a , q t 3 , • • •, q/ T , p, ..., 7 e A ( 60 ) 

where q a is the displacement vector of DNS node a, and A is the set of all DNS nodes. 
Likewise, the micromorphic finite element nodal displacements and micro-displacements are 
written as 
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particulate 

material 

target 


DNS (DE/FE/CSE) 


coupling 
between DNS 
and 

micromorphic 
continuum FE 


micromorphic 
continuum FE 


shear/compression loading 


b dns 

Ldns-fe 

K<. _ ? QgDNS 

\b fe 


Figure 4. Two-dimensional illustration of Coupling Approach 2 between grain-FE DNS and 
micromorphic FE continuum regions. 


D = [d a ,d b ,...,d c ,<f> d ,<f) e ,...,<t> f ] T (61) 

a,b ,..., c G J\f , d,e ,..., / E M 

where d a is the displacement vector of node a, <p d is the micro-displacement tensor in vector 
form of node d, A f is the set of all nodes, and A4 is the set of micromorphic finite element 
nodes with micro-displacement tensor degrees of freedom, where A4 C A f. 

3.2.2 3D Kinetic and Potential Energy Partitioning and Coupling 

For the DNS equations, the kinetic energy is T®, dissipation function F®, and potential 
energy U Q , such that 


T q = i QM q Q 

F q = a Q T Q (62) 

rQ 

U q (Q) = / F 1NT ’ Q (S)dS 

Jo 

where is the mass matrix, and F 1NT is an internal force vector for the DNS region. 
For the micromorphic continuum FE equations, T D is the kinetic energy, F D the dissipation 
function, and U D the potential energy, such that 
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T d = -DM d D 
2 

F° = cl d T d (63) 

f D 

U d (D) = / F INT ’ D (S)dS 

Jo 

where M J> is the mass and micro-inertia matrix, and p IN1 - D i s the internal force vec¬ 
tor. We assume the total kinetic and potential energy and dissipation of the coupled 
DNS-micromorphic-continuum system may be written as the sum of the energies using the 
scalar distance function u(X ) = f(X)/£ 0 , where f(X) is the orthogonal distance of point 
X E J3 dns ~ fe from dB FE , and f 0 is the orthogonal length of Jg DNS ~ FE ^ such that 


t(q,d) 

= (1 -lu)T q (Q)+ujT d (D) 

(64) 

U(Q,D) 

= (1 -u)U q (Q)+uU d (D) 

(65) 

F(Q,D) 

= (1-c o)F q (Q)+uF d {D) 

(66) 


Furthermore, the potential energy U of the coupled system is modified to include a constraint 
using a Lagrange multiplier, such that the Augmented-Lagrange multiplier constrained po¬ 
tential energy U\ ug Lag is written as 

U A ugLag = u + A • g + -kg ■ g (67) 

where A is the vector of Lagrange multipliers (forces acting on DNS nodes in the overlap 
region), g is the vector of constraints applied at each DNS node a within the overlap region, 
and k is a penalty parameter. The constraint at DNS node a in the overlap region is written 
as 

9a = ( K) h (t ) - 9a(t) = 0 ( 68 ) 

where (u' a ) h (t) is the micro-continuum displacement at node a shown in Eq.(14). Lagrange’s 
equations may then be stated as 

d_ / &T\ _ dT_ dF/ dT/AugLag _ _ jpEXT,Q 

dt \dQ) dQ OQ dQ 

d_ fdT \ _ &T_ dF/ <9P A u g Lag _ zpEXT,D 

Jt\-5D)~'dD + W + ~dD^^ F m 

which lead to a coupled system of governing equations (linear and first moment of momen¬ 
tum) for the coupled DNS-micromorphic-continuum mechanics. After taking the derivatives, 
and assuming mass and micro-inertia lumped matrices (for explicit dynamics), we arrive at 
the coupled system of equations to solve at DNS node a 

A {(1 - “(VJ) + S‘F q (Q)} -K- kg a = /f x ™} (70) 

a 

where m® is the lumped mass at DNS node a, f I a NT,( ^ is the internal force vector at DNS 
node a, f^ XT,( ^ is the external force vector at DNS node a. For micromorphic FE node a, 
we have 
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A 


m°d a 


+ a D m f d a + f 


INT,D, 
a 


D) 


+ a < A *+*9j] = t 


EXT,D 
a 


(71) 

where is the lumped mass at micromorphic FE node a, f J a NT ’ D is the internal force vector 
at micromorphic FE node a, fa XT ’ D is the external force vector at micromorphic FE node 
a. For micromorphic FE node d, we have 


LU 


:x n 


id 0 d + a Did 4>d + f 


■D . 


f INT,D, 
d ' 


D) 


A [Nd (X Q ) (A a + kg a ) <g) 3„] = 


nEXT,D 

J d 


(72) 


where is the lumped micro-inertia at micromorphic FE node d. 

An explicit time integration procedure and solution method is described in Xiao and 
Bclytschko [2004], which can also be generalized for implicit time integration and a Newton- 
Raphson nonlinear solution procedure. 
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4 Overview of adaptivity approach 

Numerical adaptivity is introduced to track, at the grain-scale, regions of dynamic material 
failure. An adaptive scheme can be considered in two ways to calculate conversion criteria: 
(a) using a calculated shear/compaction strain energy density averaged over a RVE of DNS 
micro-structure, and calculated at a Gauss point in a micromorphic FE; and (b) a posteriori 
model error estimators (e.g., Oden et al. [2006]). In terms of method development for 
conversion, we consider the cases illustrated in Fig.5 where a continuum micromorphic FE 
mesh region is converted to a particulate micro-structural DNS region, and vice versa. One 
of the challenges with this conversion is that we do not know the material micro-structure 
to make up the DNS region, nor its material parameters. We can use probability density 
functions (PDFs) to sample approximate grain size, shape, percentage of matrix versus grain 
material, inter-/intra-granular constitutive parameters, and grain-matrix interface de-bond 
strength, for instance. The sampling of PDFs will provide an initial guess for the DNS 
micro-structural geometry and material parameters. To ensure smooth conversion in terms 
of average stress, strain, and energy, we can calibrate inversely the material parameters (using 
the PDF sampling as an initial guess, and the DNS geometry fixed) through a micromorphic 
homogenization procedure. 




Figure 5. Adaptivity: conversion from micromorphic continuum FE to particulate DNS material 
model representation, and vice versa. The micromorphic continuum FE region could be the blue- 
dashed box at the bottom-center of the illustration in Fig.f to be converted to a particulate DNS 
region for adaptive high-fidelity material modeling as the projectile penetrates the target. Note that 
since we do not know the exact particulate DNS configuration and material properties when con¬ 
verting from continuum FE to DNS in a region of potential fracture nucleation and fragmentation, 
this process is probabilistic, and thus PDFs will be used for the conversion. 

We ended up not having time to implement these conversions illustrated in Fig.5. 
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5 Future Research 

Future research entails the following: (1) extend the micromorphic 3D FE implementation 
to explicit dynamics in Tahoe; (2) implement the overlap procedures presented in Sect.3; 
and (3) test the adaptivity concepts presented in Sect.4. 

Overall, the full 3D finite strain micromorphic coupled FE implementation and Drucker- 
Prager model formulation and implementation took longer than expected, but now it is 
nearly completed for the quasi-static case. The overlap coupling in 3D will continue after 
this project, to be supported by another project and to include the multiphase (solid, liquid, 
gas) nature of geologic and particulate materials. The adaptivity of our multiscale modeling 
approach will be investigated as part of future work, and it is still an important problem to 
investigate. 
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